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1.0  INTRODUCTION 


Ionospheric  mapping  and  modeling  relies  on  measured  ionospheric 
characteristics.  The  U.S.  Air  Force  DISS  (Digital  Ionospheric  Sounding 
Systems)  network  of  ionosondes  provides  real  time  electron  density  profiles 
for  the  locations  of  the  sounders.  Valuable  additional  information  can  be 
obtained  for  the  ionospheric  mid-points  between  two  stations  separated  by 
1000  to  3000  km  by  the  use  of  oblique  ionograms.  Since  the  Digisonde  256 
sounders  used  in  the  DISS  network  are  capable  to  operate  in  the  oblique 

ionogram  mode,  the  Phillips  Laboratory  initiated  the  DORIS  (Digital  Oblique 

Remote  Ionospheric  Sensing)  program  in  1987  to  establish  the  feasibility  of 
such  an  approach. 

In  a  number  of  reports,  the  University  of  Massachusetts  Lowell,  Center  for 
Atmospheric  Research  (ULCAR)  documented  the  system  requirements,  test 
data,  and  results  of  oblique  sounding  between  existing  DISS  modes 

(Kuklinski  et  al.,  1987).  Operational  implementation  of  the  DORIS  technique 
relies  on  the  development  of  software  that  automatically  inverts  the  oblique 
ionograms  into  mid-point  electron  density  profiles.  The  inversion  from 
oblique  ionogram  to  mid-point  profile  (DORIS  inversion)  consist  of  two  major 
sections  of  about  equal  difficulty:  one.  the  autoscaling  and  identification  of 
the  various  echo  traces,  and  two,  the  inversion  of  the  trace  into  a 

representative  mid-point  electron  density  profile.  The  autoscaling  section 
was  described  by  Kuklinski  et  al.,  (1992).  The  present  report  describes  an 
innovative  inversion  technique. 

During  the  period  30  April  1990  through  31  September  1992,  L’LCAR 
implemented  and  tested  a  complete  oblique  ionogram  autoscaling  and 
inversion  algorithm  using  both  synthetic  and  actual  oblique  ionogram  data. 
The  DORIS  program  uses  a  three  layer  quasiparabolic  model  to  represent  the 
midpoint  electron  density  profile.  The  inversion  process  was  formulated  as 
a  two  stage  configurational  optimization  problem,  that  was  solved  using  a 
simulated  annealing  optimization  technique.  Results  of  this  work  were 
reported  in  a  paper  (Reinisch  et  al.,  1991)  and  presented  at  the  AGARD 
Conference  in  Turkey  in  October  1991.  This  paper  is  attached  in  Appendix 
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A  to  this  report.  Section  2  reviews  some  of  the  work  reported  earlier 

(Kuklinski  et  al.,  1988,  1989,  1990). 

2.0  REVIEW  OF  OBLIQUE  lONOGRAM  NOISE  SUPPRESSION  AND 
AUTOSCALING  ALGORITHM 

The  automatic  oblique  ionogram  noise  suppression  and  scaling  algorithms 
use  a  scaled  vertical  ionogram,  produced  by  the  Digisonde  ARTIST  (Reinisch, 
et  al.,  1989),  to  synthesize  the  echo  traces  that  would  be  measured  over  the 
oblique  link  if  the  midpoint  ionosphere  and  the  ionosphere  at  the 

transmitter  were  identical.  These  synthesized  oblique  echo  traces  are  used 
by  the  subsequent  noise  suppression  filter  and  the  trace  identification 
algorithm. 

The  first  operation  performed  on  the  received  oblique  ionogram  is  a 

frequency  redundancy  noise  suppression.  A  100  kHz  oblique  ionogram  is 
generated  from  either  25  or  50  kHz  sounding  data  by  selecting  the  highest 

signal-to-noise  ratio  channel  from  each  set  of  four  adjacent  25  kHz 

soundings  or  two  adjacent  50  kHz  soundings.  The  effectiveness  of  this 
technique  is  based  on  the  fact  that  the  interfering  signals  produced  by 
spurious  radio  transmitters  are  manifested  as  large  amplitude  signals  at 
every  group  path  for  a  narrow  frequency  band;  generally  an  interfering 
transmitter  appears  at  only  one  sounding  frequency.  Since  the  group  path  is 
essentially  constant,  that  is  it  changes  less  than  the  resolution  of  the 
measurement  technique  (5  km),  for  several  25  or  50  kHz  sounding 
frequencies,  the  individual  sounding  frequency  with  the  highest  signal  to 
noise  ratio  within  a  predetermined  frequency  range  can  be  used  to 
represent  that  entire  range. 

The  subsequent  noise  suppression  filter  uses  a  minimum-mean-square 
Wiener  filtering  technique.  The  preproccssed  100  kHz  ionogram  is  modelled 
as  the  sum  of  a  noise/interference  signal  and  a  desired  noise  free  oblique 
ionogram.  The  filter  is  obtained  by  dividing  the  power  spectral  density  of 
the  synthesized  oblique  echo  traces  by  the  sum  of  the  power  spectral 
densities  of  the  synthesized  oblique  echo  traces  and  an  estimate  of  the  noise. 
The  noise  estimate  is  calculated  from  the  average  of  a  large  number  of 
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ionograms  obtained  with  the  receiver  subsystem  operating,  but  with  the 
transmitter  off.  The  2-D  Fourier  transform  of  the  preprocessed  100  kHz 
ionogram  is  multiplied  by  the  resulting  Wiener  filter.  The  inverse  2-D 
Fourier  transform  of  the  resulting  product  produces  the  potential  oblique 
echo  trace  segments.  This  process  is  outlined  in  l-igure  I. 

The  synthesized  oblique  echo  traces,  derived  from  the  ARTIST  scaled 
vertical  data,  are  used  to  define  regions  of  interest  for  subsequent 
processing  of  the  potential  oblique  echo  traces.  These  regions  of  interest  are 
defined  as  all  points  within  a  fixed  number  of  range  bins  and  sampling 
frequencies  of  the  synthesized  oblique  echo  traces.  At  present  the  allowable 
variations  were  derived  from  the  variances  of  an  oblique  ionogram  data 
base.  Any  data  in  the  region  of  interest  for  a  particular  layer  is  considered 
to  be  a  potential  part  of  the  echo  trace  for  that  layer. 

The  echo  trace  segments  in  a  particular  region  of  interest  are  scaled  using  a 
tracking  algorithm.  This  algorithm  takes  and  finds  the  best  fit  of  a  template 
of  the  echo  trace  cross-section  to  an  initial  point  on  the  echo  trace  segment. 
To  obtain  the  best  fit,  the  template  is  translated  and  rotated  about  the  initial 
point.  The  results  of  the  fitting  process  are  the  locations  of  the  center, 
leading  and  trailing  edges  of  the  echo  trace.  The  algorithm  uses  these  points 
to  determine  the  next  point  on  the  echo  trace  where  the  template  fitting 
procedure  is  repeated.  The  entire  echo  trace  is  tracked  in  this  manner. 

The  leading  edge  points  determined  by  the  tracking  algorithm  are  used  to 
determine  the  parameters  of  a  three  layer  quasiparabolic  electron  density 
profile  model. 
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Oblique  lonoaram:  s(t,f)  =  n(t,f)  +  et(t,f) 


n(t,f) 


Assumptions:  n(t,f)  Stationary  Random  Process  of 

Known  PSD  =  Pp(a)t,(Df) 

et{t,f)  Stationary  Random  Process  of 
Known  PSD  =  Pgj((0|,co^) 


Pet 


A 

Result:  et(t,f)  MMSE  Estimate  of  et(t,f) 

Figure  1.  Two-Dimensional  Wiener  Filtering 
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3.0  MODIFICATIONS  AND  IMPROVEMENTS  MADE  TO  INVERSION 

ALGORITHM 


The  most  significant  modification  of  the  previously  proposed  inversion 
algorithm  is  the  utilization  of  a  simulated  annealing  optimization  technique 
to  perform  the  inversion  process.  A  major  limitation  of  the  previously 
proposed  inversion  process,  in  particular  for  long  ground  ranges  and 

multiple  layer  data,  was  the  numerical  stability  of  the  iterative  Newton- 
Raphson  technique.  An  additional  limitation  of  this  technique  was  that 

while,  in  principle,  each  layer  could  be  described  by  more  than  one 

quasiparabolic  segment,  no  straight  forward  technique  was  available  to 
determine  the  optimum  number  of  segments  for  any  given  scaled  echo 

trace.  Typically  the  optimum  will  be  a  compromise  between  a  single 
quasiparabolic  layer,  which  may  not  adequately  describe  the  actual  layer, 

and  a  large  number  of  quasiparabolic  segments,  which  yield  a  numerically 
unstable  inversion  problem.  The  formulation  of  the  inversion  problem  as  a 
configurational  optimization  problem  allowed  a  simulated  annealing 
technique  that  yielded  the  desired  optimal  representation  of  each  layer,  to 
be  used  (Dyson  and  Bennett,  1988). 

An  additional  change  that  was  made  to  the  inversion  algorithm  was  the 
inclusion  of  a  two-stage  inversion  algorithm  that  iteratively  adjusted  the 

quasiparabolic  layer  parameters  to  satisfy  the  observed  propagation  mode 
extremes  before  the  final  simulated  annealing  optimization  technique  was 

applied.  The  two-stage  inversion  algorithm  was  implemented  by  the 

addition  of  an  iterative  procedure,  outlined  in  Figure  2.  This  procedure 
was  applied  in  sequence  to  the  E,  FI  and  F2  layers.  For  each  la>er  an 
estimate  of  the  QP  layer  parameters  were  obtained  from  either  a  scaled 
vertical  end-point  ionogram  and/or  an  ionospheric  model.  These  initial 
estimates  arc  labeled  y”,  .  r'J,  and  f”.  in  the  upper  panel  of 
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Figure  2  which  shows  the  iterative  procedure  for  a  single  layer.  In  Figure  2. 
fl  is  the  MUF  frequency,  f2  the  lowest  frequency  where  a  high  ray 
propagation  mode  was  observed  in  the  oblique  ionogram  and  fi  the  lowest 
frequency  where  a  low  ray  propagation  mode  was  observed  in  (he  oblique 
ionogram.  The  procedure  adjusted  the  critical  frequency  of  the  layer  until 
no  ray  reached  the  known  ground  range  for  any  frequencies  above  the 
observed  MUF  frequency,  as  seen  in  the  upper  panel  of  Figure  2;  adjusted 
the  semithickness  of  the  layer  until  at  least  one  high-angle  ray  reached  the 
known  ground  range  at  the  lowest  observed  high-angle  ray  sounding 
frequency  as  seen  in  the  middle  panel  of  Figure  2;  and  adjusted  the  base 
height  of  the  layer  until  at  least  one  low-angle  ray  reached  the  known 
ground  range  at  the  lowest  observed  low-angle  ray  sounding  frequency,  as 
seen  in  the  lower  panel  of  Figure  2.  This  sequence  of  adjustments  was 
repeated  until  all  three  of  the  observed  propagation  mode  extreines  were 
simultaneously  satisfied.  This  initial  "coarse"  inversion  process  was  tested 
using  the  Wallops  Island  to  Goose  Bay  oblique  ionograms.  In  every  case  the 
procedure  rapidly  converged  to  a  quasiparabolic  layer  model  that  produced 
group  paths  at  all  sounding  frequencies  that  were  within  50  km  of  the 
actual  group  paths.  The  subsequent  simulated  annealing  optimization 
algorithm  made  "fine"  adjustments  in  the  layer  parameters  resulting  in 
electron  density  profiles  with  critical  frequencies  within  200  kHz  of  the 
correct  values  and  peak  heights  within  20  km  of  the  correct  values. 

The  simulated  annealing  optimization  technique  used  an  analogy  between 
the  configurational  optimization  problem  of  interest  and  the  cooling  of  a 
thermodynamic  system.  The  energy  of  a  thermodynamic  system  that  has 
been  slowly  cooled  is  minimum,  while  a  quickly  cooled  system  typically 
arrives  at  a  higher  final  energy  state.  These  are  analogous  to  the  desired 
global  minimum  and  a  local  minimum,  respectively,  of  the  given 
configurational  optimization  problem.  The  simulated  annealing  technique 
avoids  the  local  minima  produced  by  most  iterative  optimization  techniques 
that  employ  some  type  of  gradient  algorithm  by  allowing  changes  in 
configuration  that  may  increase  the  desired  cost  function  for  some  fraction 
of  the  iterations.  The  cost  function  that  was  minimized  in  the  electron 
density  inversion  problem  was  the  total  squared  error  between  the 
measured  group  path  vs.  sounding  frequency  data  and  the  group  path  vs. 
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sounding  frequency  data  produced  by  the  quasiparabolic  model.  The 
configuration  of  the  model  consisted  of  the  specific  values  of  semithickness, 
peak  height  and  critical  frequency  for  each  layer. 

4.0  RESULTS 

The  completed  inversion  algorithm  was  first  tested  using  synthetic  data 
generated  by  the  three  layer  quasiparabolic  model  analytic  ray  tracing 
software  developed  during  this  project.  A  typical  synthetic  oblique 
ionogram  is  seen  in  Figure  3.  This  ionogram  was  calculated  for  a  ground 
range  of  2060  km,  the  length  of  the  Goose  Bay  to  Wallops  Island  link  used 
to  produce  actual  oblique  ionograms.  The  analytic  ray  tracing  software 
determined  both  one  and  two  hop  propagation  modes  and  identified  echo 
trace  points  as  either  high  or  low  ray  modes.  This  identification  is 
important  during  the  simulated  annealing  optimization  process.  Since  a 
complete  QP  electron  density  profile  model  for  each  layer  allows 
propagation  modes,  for  a  given  sounding  frequency  and  ground  range,  an 
ambiguity  can  exist  during  the  inversion  process.  By  tugging  each  point  on 
the  oblique  ionogram,  the  ambiguity  is  removed  and  the  inversion  process 
can  be  more  efficient  and  accurate. 

Using  the  echo  trace  points  indicated  with  circular  dots  in  Figure  3,  the 
inversion  algorithm  produced  the  profile  seen  in  Figure  4.  The  electron 
density  used  to  generate  the  oblique  ionogram  and  the  electron  density 
produced  by  the  inversion  algorithm  are  essentially  identical. 

Figures  and  6  illustrate  the  application  of  the  DORIS  algorithms  to  real 
data.  Figure  .'Sa  shows  a  typical  Goose  Bay  to  Wallops  Island  oblique 
ionogram.  The  ±  symbols  mark  the  autoscaled  trace  points.  Figure  6 
shows  the  midpoint  electron  density  profile  produced  by  the  DORIS 
inversion  algorithm  together  with  the  mid-point  electron  density  profile 
determined  from  a  simultaneous  vertical  ionogram  recorded  at  the  Maine 
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Figure  4.  Inversion  Using  Data  from  Figure  3 
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Figure  6.  Inversion  Using  Data  from  Figure  5 
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OTH  radar  site.  The  electron  density  profile  produced  by  the  DORIS 
inversion  matches  the  actual  electron  density  profile  quite  well  near  and 
below  the  peaks  of  the  FI  and  F'2  layers.  The  deep  "valleys"  in  the  DORIS 
profile  result  from  the  present  limitation  of  one  quasiparabolic  segment 
per  layer.  The  difference  in  the  critical  frequencies  between  the  measured 
midpoint  eleetron  density  profile  and  the  DORIS  profile  occurred  because 
no  high  ray  E  layer  points  were  detected  by  the  scaling  algorithm.  All 
scaled  E  layer  points  were  from  the  region  10-20  km  below  the  peak 
height.  Since  no  constraints  were  placed  on  the  E  layer  critical  frequency, 
the  optimization  algorithm  found  that  the  best  match  to  the  observed  low 
ray  E  layer  group  path  points  was  obtained  using  a  quasiparabolic  layer 
with  a  critical  frequency  200  KMz  greater  than  the  actual  E  layer  critical 
frequency.  The  inclusion  of  a  multisegment  representation  of  each  layer  in 
an  advanced  version  of  DORIS  would  have  corrected  this  type  of  error. 

3.0  STATUS  OF  EXISTING  PROCESSING  SOFTWARE 

The  present  DORIS  software  package  that  produces  midpoint  electron 
density  profiles  from  oblique  ionograms  is  written  in  Fortran.  The  noise 
suppression  and  trace  enhancement  Wiener  Filter  algorithm  and  the  echo 
trace  identification  algorithm  was  implemented  on  an  Apollo  3500  work 
station.  This  algorithm  requires  a  scaled  end-point  vertical  ionogram  and 
noise  power  spectral  density  to  process  an  oblique  ionogram.  The  filtering 
and  trace  identification  takes  approximately  50  seconds  per  oblique 
ionogram  on  the  Apollo  platform.  The  inversion  algorithm,  also  written  in 
Fortran,  was  implemented  in  a  DOS  environment  on  a  PC  386.  The 
inversion  algorithm  consists  of  a  main  program  and  two  subroutines:  the 
first  is  the  three  layer  quasiparabolic  ray  tracing  module  and  the  second 
computes  the  takeoff  angle  needed  to  satisfy  a  given  group  path  and 
sounding  frequency  point  for  a  specified  ground  range.  The  source  code 
for  these  algorithms  is  available  on  floppy  disk.  Appendix  B  contains  an 
operator’s  manual  for  the  DORIS  inversion  program.  The  inversion 
software  requires  absolute  group  path  vs.  sounding  frequency  data  labeled 
as  to  which  echo  trace  they  are  from  and  whether  they  are  high  or  low  ray 
propagation  modes,  and  a  three  layer  quasiparabolic  estimate  of  the 
midpoint  electron  density  profile.  On  a  386-based  PC  system  the  current 
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version  of  the  inversion  algorithm  takes  between  2  and  8  minutes  to 

determine  the  optimal  midpoint  electron  density  profile.  We  estimate  a 
total  DORIS  run  time  on  a  386  PC  of  between  4  to  10  minutes,  when  the 
autoscaling  and  profile  inversion  algorithms  are  combined. 

6.0  CONCLUSION 

The  results  of  the  activities  conducted  during  the  DORIS  project  have 
produced  a  complete  prototype  software  package  for  automatic  scaling  and 
inversion  of  oblique  ionograms.  The  system  has  been  designed  to  utilize 
any  available  information  about  the  desired  midpoint  electron  density 
profile.  Quantitative  information  derived  from  endpoint  vertical 
ionograms  and/or  ionospheric  models  is  used  in  the  scaling  and  inversion. 

The  noise  cleaning  and  echo  trace  enhancement  Wiener  filtering  method 
produced  good  estimates  of  the  echo  traces  present  in  ail  of  the  recorded 
oblique  ionograms.  The  conversion  of  the  present  off-line  version  of  this 
algorithm  to  an  on-line  PC  based  implementation  would  be  the  only 

modification  necessary  in  an  operational  system. 

The  simulated  annealing  optimization  algorithm  produces  good  estimates 

of  the  actual  midpoint  electron  density  profiles  if  at  least  3  to  5  points  on 

both  the  high  and  low  rays  of  each  echo  trace  are  available.  The  time 

required  for  the  present  inversion  algorithm  to  determine  a  midpoint 

electron  density  profile  could  be  reduced  if  the  annealing  schedule  portion 
of  the  optimization  algorithm  was  improved.  If  the  control  parameter 
could  be  reduced  at  an  optimal  rate  such  that  the  algorithm  would  be 
better  able  to  recognize  local  minima  during  the  optimization  process  and 
produce  the  midpoint  electron  density  profile  faster.  The  general  question 
of  how  to  determine  an  optimal  annealing  schedule  is  still  an  active 
research  topic  and  hence,  the  optimal  annealing  schedule  for  the  oblique 

inversion  optimization  method  will  require  additional  effort. 

The  present  echo  trace  scaling  algorithm  is  the  portion  of  the  total  DORIS 
processing  program  that  limits  overall  performance.  The  regions  of 
interest  around  each  echo  trace  that  limit  the  extent  of  the  search 
algorithm  often  exclude  portions  of  the  echo  trace  of  interest  or  they 


contain  more  than  one  echo  trace.  A  preliminary  study  of  an  alternative 
method  of  trace  identification  based  on  an  attributed  relational  graph 
(ARG)  representation  of  oblique  ionograms  has  been  conducted.  This 
method,  which  has  been  proposed  for  other  machine  vision  tasks, 
represents  an  oblique  ionogram  as  a  graph  (Eshera  and  Fu,  1984).  The 
nodes  of  the  graph  represent  the  primitive  elements  of  the  oblique 
ionogram,  which  in  our  work  are  a  set  of  straight  line  segments.  The 
branches  of  the  graph  represent  the  spatial  relationships  between  the 
primitive  elements.  The  oblique  ionogram  echo  traces  are  found  by  first 
filtering  the  raw  oblique  ionogram  to  enhance  the  region’s  that  contain 
potential  echo  trace  straight  line  segments.  The  ARG  representations  of  all 
possible  combinations  of  these  potential  echo  trace  segments  are  compared 
with  a  reference  ARG  obtained  from  the  endpoint  electron  density  profile. 
The  echo  trace  for  each  layer  is  produced  from  the  subset  of  all  potential 
echo  trace  segments  that  best  matches  the  reference  graph.  The  method 
allows  both  the  general  structure  and  position  of  each  echo  trace  and  the 
specific  data,  such  as  expected  critical  frequencies,  to  be  used  in  trace 
identification  process.  The  refinement  and  implementation  of  this  trace 
identification  method  would  produce  a  more  robust  automated  oblique 
ionogram  scaling  and  profile  inversion  algorithm. 
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SUMMARY 

.A  small  low  cost  digital  ionosonde,  the  Digisonde 
Portable  Sounder  (DPS)  has  been  developed 
which  uses  500ps,  10%  duty  factor,  wide  pulses 
for  vertical  sounding  and  8.5ms  pulses  for 
oblique  sounding.  Intrapulse  coding  and  pulse 
compression  techniques  result  in  a  67ns 
resolution  for  both  waveforms.  A  new 
autoscaling  technique  for  oblique  ionograms 
inverts  the  oblique  echo  traces  into  midpoint 
electron  density  profiles  that  arc  modeled  as  a 
sum  of  quasiparabolic  layers. 

1 .  INTRODUCTION 

Reliable  high  frequency  (HE)  communication. 
Over-the-Horizon  radar  tracking,  and  HF  direc¬ 
tion  finding  all  depend  on  good  real  lime 
knowledge  of  the  ionosphere.  In  the  last  decade 
modern  digital  ionosondes  (Reinisch. 1986)  be¬ 
came  available  that  automatically  analyze  the 
vertical  ionograms  and  provide,  in  real  time,  the 
vertical  electron  density  profiles  for  the  location 
of  the  sounder  (Reinisch  et  al.,  1990).  For  many 
HF  applications  the  ionosphere  must  be  specified 
in  a  large  area  and  a  network  of  ionosondes  is 
required.  The  new  Digisonde  Portable  Sounder 
(DPS)  has  been  designed  for  optimal  efficiency  in 
such  a  network  by  providing  the  necessary- 
hardware  and  software  for  both  vertical  and 
oblique  (bistatic  between  DPS  stations)  sounding. 

For  vertical  sounding  the  transmitter  pulse 
width  is  limited  to  about  500txs  considering  the 
pulse  propagation  time  of  600ns  for  E  region 
echoes.  Intrapulse  phase  coding  and  pulse  com¬ 
pression  techniques  provide  a  range  resolution 
of  67ns  (or  10  km).  For  quasi-vertical  or  bistatic 
sounding  the  DPS  transmits  8.5ms  pulses  at 
100%)  duly  cycle  with  the  same  range  resolution 
of  67ns.  A  DPS  network  can  accurately  map  the 
ionosphere  if  both  vertical  and  oblique 
ionograms  are  used  to  specify  the  electron 
density  profiles  at  the  locations  of  the  sounders 
as  well  as  at  their  mid-points.  Any  successful  HF 


application  requires  "now -casting’  lor  the 
ionospheric  propagation  medium,  i.  e.  current 
specification  of  the  ionospheric  region  used  by 
the  HF  link.  Automatic  scaling  of  the  \erncal 
and  oblique  ionograms  in  real  time  provides  the 
inputs  for  the  now-casting. 

Autoscaling  of  the  vertical  ionograms  is  achieved 
with  the  ARTIST  software  (Reinisch  and  Huang. 
1983,  Gamache  et  al..  1985)  which  is  routinely 
used  in  some  forty  Digisonde  256  sounders 
worldwide  (Reinisch  et  al..  1990),  New  autoscal¬ 
ing  for  oblique  ionograms  and  inversion  to  mid¬ 
point  quasiparabolic  electron  density  profiles  is 
described  in  Section  3  of  this  paper.  Section  2 
gives  a  brief  description  of  the  Digisonde 
Portable  Sounder. 

2  .  DESCRIPTION  OF  THE  DIGISONDE 
PORTABLE  SOUNDER 

Ground  based  HF  observations  of  the  ionosphere 
rely  on  the  accurate  measurement  of  the  follow¬ 
ing  observables: 

1.  frequency 

2.  height/range  (time  delay) 

3.  amplitude 

4.  phase 

5.  Doppler  shift  or  spread  (frequency  disper¬ 
sion) 

6.  wave  polarization  (left  or  right  hand) 

7.  angle  of  arrival 

All  of  these  parameters  are  measured  simulta¬ 
neously  by  the  DPS  and  as  many  as  five  arc  dis¬ 
played  at  once  in  the  various  display  modes 
available.  For  instance,  the  Doppler  lonogram  in 
Figure  1  (Millstone  Hill,  4  November  199]  i 
represents  an  echo  at  a  given  fret|uenrs  and 
height  by  its  position  in  the  display.  Ihe 
a  m  n  I  i  1  u  d  e  (in  2dB)  is  represented  numerically 
by  an  optically  weighted  font,  the  Doppler  shift 
is  represented  by  shades  of  color  while  the 
polarization  is  represented  by  the  color  group. 
This  paper  shows  a  black/white  prinioui  using 
the  letter  X  to  indicate  X-echoes  Ihe 
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F'lgure  1.  Dictsonde  Pori;iblc  Sounder  Vertical 

ionogram;  Millstone  Hill.  Massachusetts 

auiomaticail y  scaled  leading  edges  of  the  E  and 
F-traces  are  identified  by  the  letters  E  and  F,  and 
the  electron  density  profile  is  indicated  by  the 
letter  T. 

Accurate  measurement  of  all  seven  parameters, 
except  the  first,  depends  critically  on  the  signal- 
to-noise  ratio  (SNR)  of  the  received  signal, 
therefore  vertical  incidence  ionospheric  sounders 
capable  of  acquiring  high  quality  scientific  data 
have  hisioncally  utilised  powerful  pulse 
transmitters  in  the  2  to  30l:Watt  range.  For 
instance,  measuring  phase  to  I  degree  accuracy 
requires  a  SNR  better  than  40dB  (assuming  a 


Gaussian  noise  distribution  which  is  actually  a 
best  case),  measuring  Doppler  spectra  to  ICi 
amplitude  accuracy  requires  over  20dB  SNR,  and 
precision  ranging  by  the  do/df  technique 
requires  better  than  20dB  SNR.  Making  a  small 
inexpensive  portable  sounder  requires  use  of  a 
small  low  power  transmitter  which  makes  it 
difficult  to  achieve  the  required  SNR.  However, 
advanced  digital  signal  processing  techniques  in 
the  DPS  overcome  this  difficulty. 

Physical  Layout  of  the  DPS 

The  DPS  (Fig,  2)  is  a  software  implementation  of 
the  well  proven  measurement  techniques  used 
by  the  Digisonde  2.‘i6  sounder  (Reinisch  et  al. 
1989)  developed  some  ten  years  ago  at  the 
University  of  Lowell  Center  for  Atmospheric 
Research  (ULCAR).  Pulse  compression  tech¬ 
niques  allow  reduction  of  the  transmitter  peak 
pulse  power  from  lOkWatt  (for  the  Digisonde 
256)  to  SOOWati. 

Figure  3  shows  the  dimensions  and  component 
layout  of  the  DPS.  The  single  19"  rack  mountable 
chassis  is  installed  in  a  lightweight  fiberglass 
transport  case  which  includes  a  shock  mounted 
19"  rack.  This  unit  by  itself  (with  the  addition  of 
a  24  volt  battery  and  antennas)  is  ready  to 
make,  scale  and  store  automated  measurements 
as  soon  as  power  is  turned  on.  A  typical  instal¬ 
lation  would  also  include  a  keyboard  and  moni¬ 
tor  and  two  printers  for  real-time  color  or  black 
and  while  outputs  of  measurement  data.  The 
operator  can  control  the  system  either  locally  or 
remotely  (via  a  telephone  modem). 

Signal  Processing 

As  mentioned  above,  the  requirement  for  a  good 
SNR  seems  incompatible  with  a  small, 
lightweight,  low-cost  portable  system,  however 
by  lengthening  the  transmitted  pulse  of  a  small 
low  voltage  solid  slate  transmitter  we  can 
transmit  an  amount  of  energy  per  pulse  equal  to 
that  transmitted  by  a  high  power  pulse 
iransmiiicr  without  having  lo  provide 
components  to  handle  the  large  voltages.  The 
range  resolution  is  recovered  by  phase 
modulation  using  an  optimized  phase  code.  This 
inirapulse  modulation  allows  pulse  compression 
in  the  receiver  to  restore  a  10km  resolution  to 
separate  the  echoes  from  different  ionospheric 
layers.  The  compressed  pulse  must,  to  an 
acceptable  degree.  reproduce  a  simple 
rectangular  pulse  (i.e.  the  spurious  output  of  the 
pulse  compression  process  must  not  obscure 
other  signals  of  interest)  and  for  monostatic 
operation  (colocated  iransmiiicr  and  receiver) 
ihe  entire  wavelorm  musi  be  iransnmied  before 
ilic  first  echo  of  interest  is  received.  Therefore, 
for  vertical  incidence  sounding  we  selected  an  8 
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bn  coniplenifnur\  (Golay,  1961)  phase  code  pair 
modulated  onto  a  533usec  pulse  which  provides 
I2dB  of  SNR  improvement  upon  pulse 
compression. 

For  bistatic  operation  such  as  oblique  sounding 
or  backscatter  sounding,  a  100%  duty  cycle  CW 
(continuous  wave)  waveform  can  be  used  since 
the  transmission  need  not  be  extinguished  prior 
to  sampling  data  at  the  receiver.  For  this  mode 
of  operation  we  have  selected  the  127  bit  maxi¬ 
mal  length  sequence  code  (Sarwate  and  Pursley, 
1980)  which  allows  pulse  compression  with  a 
21dB  signal-to-noise  ratio  improvement  while 
only  producing  0.9%  spurious  amplitude  re¬ 
sponse  at  ranges  (time  delays)  other  than  the 
true  range. 

Added  to  the  12  or  21dB  of  pulse  compression 
gain  is  the  SNR  gain  offered  by  the  Doppler  inte¬ 
gration  process  (also  used  in  the  Digisonde  256) 
which  for  12S  line  Doppler  spectra  is  21dB.  The 
effect  of  all  this  signal  processing  is  tremendous. 
For  instance,  compared  to  the  SOOWatt  DPS 
oblique  waveform  with  its  42dB  SNR  improve¬ 
ment,  a  single  pulse,  non-integrating  measure¬ 
ment  system  would  require  SMWatts  of 
transmitter  posvcr  to  obtain  an  equivalent  SNR. 
Even  a  somewhat  sophis'icated  system  which 
coherently  integrates  8  pulses  (any  more  than  8 
would  require  "Digisondc-like"  Doppler  process¬ 
ing  due  to  Doppler  induced  phase  changes  over 
the  integration  period),  would  require  a  IMWait 
transmitter  in  order  to  be  equivalent. 

■Multiples  in  g 

Using  any  of  the  above  waveforms,  the  Digisonde 
Portable  Sounder  allows  multiplexed  Doppler 
integration  of  up  to  64  separate  combinations  of 
frequencies,  antennas  (the  system  includes  4 
phase  matched  receiving  antennas),  and  polar¬ 
izations  fi.e.  ordinary  or  extraordinary  propaga¬ 
tion  modes,  hereafter  referred  to  as  O  and  X). 
This  multiplexing  is  implemented  by  changing 
the  system  s  operating  parameters  from  pulse- 
to-pulse.  This  essentially  simultaneous  mea¬ 
surement  at  multiple  frequencies,  antennas  and 
polarizations  is  accomplished  by  synchronizing 
the  computers  data  acquisition  and  signal  pro¬ 
cessing  sot'ware  to  the  transmitted  waveform 
through  a  f.arriware  interrupt.  The  interrupt  is 
timed  to  o^cur  lust  250pscc  prior  to  the  occur¬ 
rence  of  a  tr.insmitted  pulse,  or  coded  waveform, 
such  that  the  fre(]uency,  antenna  sclcciion,  gain, 
polarizatioi;.  or  waveform  can  be  changed  in  a 
pre-compuied  seijuencc  The  switching  of  any  or 
all  ol  these  p.irameicrs  is  accomplished  in  22(iscc 
giving  the  ^ysIe^l  ample  time  (about  20()pscc)  to 
‘c;:le  bcNf.  ii.uiNmitimg.  .Ni  the  end  of  ilie 


multiplexed  Doppler  integration  the  system 
contains  an  entire  Doppler  spectrum  for  each  re¬ 
solvable  range,  each  receive  antenna,  each  fre¬ 
quency  and  each  polarization  measured,  up  to 
8192  128-line  complex  spectra.  The  real-time 
pulse  compression  and  Doppler  processing  is 
performed  in  a  digital  signal  processor  embed¬ 
ded  within  the  system  which  is  fast  enough  to 
allow  it  to  keep  up  in  real-time  with  the  data  ac¬ 
quisition.  The  pulse-to-pulse  control  offers  the 
flexibility  to  optimize  the  multiplexing  process 
for  any  of  a  number  of  different  types  of  mea¬ 
surements.  A  few  of  these  are: 

1.  Angle  of  arrival.  Since  the  data  on  four 
spaced  receiving  antennas  is  acquired  quasi- 
simultaneously,  the  signal  phases  measured 
on  each  of  the  antennas  can  be  used  to 
compute  an  angle  of  arrival  for  each  range 
and  Doppler  component  of  the  received 
signal. 

2.  Prf  .  Group  Height.  By  observing  the 
'hange  in  phase  of  a  signal  as  a  function  of 
small  changes  in  frequency  (e.g.  lOkHz),  the 
altitude  of  a  reflecting  layer  can  be  deter¬ 
mined  with  high  accuracy.  The  limitation  is 
signal-to-noise  ratio  dependent,  but  accura¬ 
cies  of  lOO’s  of  meters  are  routinely  achiev¬ 
able.  The  time  history  of  the  precise  height 
can  be  measured  by  repeated  soundings  at  a 
fixed  frequency. 

3.  Identification  of  O  or  X  Propagation  Mode.  By 
simultaneously  integrating  all  heights  at  both 
polarizations,  the  detected  echo  is  compared 
in  the  two  separate  data  buffers  to  determine 
which  receive  polarization  produced  a  larger 
response,  therefore  even  clliptically  polarized 
signals  can  be  declared  O  or  X  even  though 
some  of  the  echo  amplitude  ''leaks"  into  the 
wrong  polarization  channel. 

4.  High  Doppler  Resolution.  Although  high 
Doppler  resolution  could  be  achieved  simply 
by  slowing  down  the  pulse  repetition  rate  or 
producing  extremely  large  Doppler  spectra 
(e.g.  4096  complex  Doppler  lines),  both  ap¬ 
proaches  increase  the  time  required  to  make 
measurements.  Since  the  ionosphere  rarely 
produces  more  than  a  few  Hertz  of  Doppler 
shift  the  Doppler  analysis  of  pulses  transmit¬ 
ted  at  a  2{X)Hz  rate  (giving  a  Doppler  range  of 
+  l00Hz)  is  usually  inefficient.  However,  by 
simultaneously  integrating  16  sounding 
frequencies,  the  Doppler  range  is  brought 
down  to  g^GHz,  which  can  provide  0.094Hz 
Doppler  resolution  rather  than  l.5Ilz.  At  the 
end  of  this  multiplexed  integration  16  Irc- 
quency  steps  have  been  measured  and  arc 
ready  for  display.  This  high  Doppler  resolu¬ 
tion  swept  frequency  lonogram  is  performed 
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16  times  faster  than  would  be  possible  with  a 
conventional  non-multiplexed  mode  of  op¬ 
eration. 

Automatic  Sequencing 

An  operator  who  is  only  casually  familiar  with 
the  inner  workings  of  a  DPS  sounder  can  easily 
compose  some  very  complex  measurement  se¬ 
quences  by  selecting  programmed  measurement 
parameters  from  an  interactive  screen.  The  pa¬ 
rameters  presented  on  the  screen  specify  fre¬ 
quency  range,  frequency  resolution,  multiplexing 
parameters,  height  resolution,  Doppler  range  and 
resolution,  data  storage  formats  and  automati¬ 
cally  schedule  unmanned  operation  Since  there 
arc  up  to  8,192  Doppler  spectra  in  memory  at 
the  end  of  a  coherent  integration,  and  typically 
less  than  1%  is  useful  data  at  the  end  of  a 
coherent  integration,  data  storage  formats  allow 
the  user  to  specify  criteria  for  reducing  the 
amount  of  data  stored; 

a.  selection  of  1  maximum  amplitude  in  each 
Doppler  spectrum  per  height  measured. 

b.  storage  oC  spectra  from  of  a  limited  number 
of  heights  around  the  leading  edge  of  a  layer. 

c.  storage  of  a  limited  number  of  Doppler  lines 
around  zero  for  each  height  measured. 

d.  or  storage  of  everything  available  (i.e.,  a 
scattering  function  which  is  all  Doppler  lines 
at  all  ranges  measured). 

The  program  and  schedule  settings  selected  by 
the  user  (either  remotely  or  locally)  during  an 
interactive  editing  session  are  stored  in  a  hard 
disk  file  and  become  the  new  default  settings  the 
next  time  the  system  is  turned  on. 

Oblique  Soundinc 

Before  addressing  the  use  of  oblique  sounding 
data  for  the  characterization  of  the  ionosphere  at 
the  midpoint,  we  will  describe  this  operating 
mode.  In  addition  to  the  CW  waveform,  the  DPS 
has  optimized  oblique  channel  measurements  by 
providing  precision  timing,  a  highly  stable  fre¬ 
quency  source  and  automatic  synchronization 
with  distant  transmitters.  All  timing  waveforms 
are  generated  from  a  counter  chain  which  can  be 
driven  by  a  Rubidium  frequency  standard.  The 
counter  chain  can  be  synchronized  to  an 
"external”  1  PPS  (pulse  per  second)  timing  refer¬ 
ence  (typically  GPS.  GOES  satellite  clock,  or 
LORAN)  and  will  maintain  synchronization  to 
Ipscc  accuracy.  The  internal  timing  signals  can 
be  offset  by  a  pre-dciermined  amount  (i.e.  a 
menu  parameter)  up  to  0.5  seconds  to  compen¬ 
sate  for  propagation  time  delays.  Therefore, 
each  transmitted  pulse  occurs  just  as  the 
receiving  system  expects  it. 


The  oblique  lonogram  data  stored  for  midpoint 
measurements  is  a  16  line  Doppler  spectrum  ai 
each  of  256  heights  (time  delays)  providing  96 
dB  of  dynamic  range  in  0.375dB  steps.  This  data 
is  stored  on  a  disk  file  and  on  a  tape  cartridge 
and  the  midpoint  profile  processing  can  proceed. 

3.  ADAPTIVE  OBLIQUE  lONOGRA.M 

SCALING  AND  PROFILE  INVERSION 

We  have  developed  an  adaptive  set  of  robust 
algorithms  to  automatically  scale  and  invert 
oblique  ionograms.  These  algorithms  utilize  one 
or  both  scaled  end  point  vertical  ionograms  to 
estimate  the  oblique  ionogram.  This  estimate, 
which  is  used  by  subsequent  filtering  and  trace 
identification  algorithms,  is  updated  whenever 
the  vertical  ionograms  change. 

Algorithm  Development 

The  automatic  oblique  ionogram  scaling  algo¬ 
rithm  uses  a  tracking  algorithm  (Sun,  1989)  to 
determine  the  leading  edge  points  of  each  echo 
trace.  These  scaled  leading  edge  points  are  used 
directly  to  calculate  a  midpoint  electron  density. 
The  direct  inversion  process,  assumes  each  iono¬ 
spheric  layer  can  be  modeled  w-ith  a  quasi¬ 
parabolic  electron  density  profile.  This  inversion 
technique  requires  the  absolute  group  path 
information  that  is  provided  by  the  synchronized 
operation  of  two  Digisondes.  The  tracking 
algorithm  uses  the  scaled  vertical  ionogram, 
produced  in  the  end  point  Digisonde,  to 
synthesize  the  echo  traces  that  would  be 
measured  over  the  oblique  link  if  the  midpoint 
ionosphere  and  the  ionosphere  at  the  transmitter 
were  identical.  While  the  overhead  and  midpoint 
ionospheres  are  typically  not  equal,  they  are 
similar  enough  to  provide  an  excellent  starting 
point  from  which  any  differences  between  the 
ionosphere  at  these  two  locations  can  be 
oeiermineo.  inese  synthesized  oblique  echo 
traces  are  also  used  in  a  noise  suppression  filter 
and  subsequent  trace  identification  algorithm. 
The  utilization  of  these  real-time  synthesized 
oblique  echo  traces  produces  an  adaptive 
algorithm  whose  performance  is  superior  to  a 
non-adaptive  approach. 

Frequency  Redundancy  Filter 

The  first  operation  performed  on  the  received 
oblique  ionogram  is  a  frequency  redundancy 
noise  suppression.  A  lOOkllz  oblique  ionogram. 
which  has  sufficient  resolution  for  calculating  the 
equivalent  midpoint  electron  density  profile,  is 
generated  from  either  25  or  50kHz  sounding 
data  by  selecting  the  highest  signal-to-noise 
ratio  channel  from  each  set  of  four  adjacent  25 
kHz  soundings  or  two  adjacent  50kHz  soundings. 
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Noi''e  Supprevsion  Wiener  Filler 

The  subsequent  noise  suppression  filter  uses  an 
adaptive  mini  mum-mean -  square  Wiener 
filtering  technique  (Andrews  and  Hunt,  1969, 
Kuklinski  et  al.,  1990).  The  lOOkHz  oblique 
ionogram  is  modeled  as  the  sum  of  a 
noise/interference  signal  and  the  desired  noise 
free  oblique  autotraces.  The  Wiener  filter  is 
produced  by  dividing  the  power  spectral  density 
of  the  synthesized  oblique  echo  traces  by  the 
sum  of  the  power  spectral  densities  of  the 
synthesized  oblique  echo  traces  and  an  estimate 
of  the  noise  power  spectral  density  as  shown  in 
equation  (It. 

Pet(  wt.tof) 

H  (  cot.  wf  )  = -  ( 1 ) 

Pei(  (Uptof)  +  Pn( 

The  noise  power  spectral  density  was  calculated 
from  ionograms  recorded  while  the  transmitter 
was  off,  using  a  two-dimensional  periodogram 
technique.  Typical  results  are  shown  in  Figure  4. 
The  upper  panel  is  a  representative  lOOkHz 
oblique  ionogram  with  the  corresponding  filtered 
ionogram  tn  the  lower  panel. 

The  synthesized  oblique  echo  traces,  derived 
from  the  end  point  ionograms,  are  also  used  to 
define  regions  of  interest  for  subsequent 
processing.  These  regions  of  interest,  as  shown 
in  the  upper  panel  of  Figure  5,  are  defined  as  all 
points  within  a  specified  number  of  group  path 
range  bins  and  sounding  frequencies  of  the 
synthesized  oblique  echo  traces.  At  present  the 
allowable  variations  were  derived  from  the 
variances  of  an  oblique  ionogram  data  base.  Any 
echoes  in  the  region  of  interest  for  a  particular 
layer  is  considered  to  be  part  of  the  echo  trace 
for  that  layer. 

The  tracking  algorithm  determines,  at  a  given 
sounding  frequency,  where  the  echo  trace  most 
closely  matches  an  assumed  template  of  echo 
trace  amplitudes.  To  obtain  this  optimal  location, 
the  template  is  translated  and  rotated  about  an 
initial  estimate.  The  mean  square  difference 
between  the  actual  echo  trace  amplitudes  and 
the  template  is  used  by  two  sequential  gradient 
searches,  cne  for  translation  and  one  for  rotation, 
to  determine  the  optimal  location.  The  results  of 
this  local  search  process  are  the  locations  of  the 
center,  leading  and  trailing  edges  of  the  echo 
trace.  The  algorithm  uses  this  information  to 
estimate  the  group  path  of  the  echo  trace  at  the 
next  sounding  frequency  to  be  processed.  The 
template  matching  procedure  is  repeated  until 
the  enure  echo  trace  is  scaled  in  this  manner. 


MiJrioint  Profile  Modeling 

Our  technique  assumes  that  the  ionosphere  at 
the  midpoint  can  be  modeled  by  one  or  more 
quasiparabolic  (QP)  segments  (Dyson  and 
Bennett,  1988).  In  principle,  the  parameters  that 
define  a  single  QP  layer  can  be  calculated  from 
any  three  pairs  of  group  path  vs.  sounding 
frequency  data  since  the  ground  range  between 
transmitter  and  receiver  is  known.  in  practice 
however,  the  numerical  sensitivity  of  the 
equations  that  relate  the  group  path  and  ground 
range  of  a  ray  to  the  take-off  angle  and  QP  layer 
parameters,  causes  these  three  point  solutions  to 
be  unstable.  The  tracking  algorithm  utilized  here 
has  the  flexibility  to  produce  the  required 
number  of  scaled  leading  points  for  a 
numerically  stable  inversion. 

The  inversion  algorithm  uses  a  simulated 
annealing  optimization  method  to  determine  the 
QP  electron  density  profile  (Croft  and  Hoogasian, 
1968)  for  each  layer.  The  QP  layer  is  defined  as  ; 


I  0  elsewhere 


where  Ne=  electron  density,  r  =  radial  distance 
from  earth  center,  rb  =  radial  distance  to  the 
bottom  of  the  layer,  rm  =  radial  distance  of 

maximum  electron  density,  =  maximum 

electron  density,  and  y,„  the  layer  semithickness. 

The  E  layer  QP  parameters  are  determined  from 
the  scaled  oblique  E  echo  points.  This  QP  E  layer 
and  the  scaled  oblique  FI  echo  points  are  used  to 
determine  the  FI  layer  QP  parameters  and,  if  an 
F2  trace  is  present,  the  E  and  FI  layer  QP  pa¬ 
rameters  and  the  scaled  oblique  F2  echo  points 
are  used  to  determine  the  F2  layer  QP 

parameters. 

The  multilayer  QP  model  produces  the  following 
analytical  relationships  for  the  group  path,  P', 
and  the  ground  range,  D,  of  a  ray  in  terms  of 
takeoff  angle,  sounding  frequency,  and  layer 
parameters. 

D=  gli(  tmi' ymi- Nmi.  P>0  (2) 

P  =  ^  g2i(  ''mi'  ymi-  Nmi,  p,f)  (3) 

where  i  is  the  layer  index,  P  is  the  ray  takeoff 
angle,  and  f  is  the  sounding  frequency. 
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Figure  4.  lOOkHz  Oblique  lunogrum  Produced  via  the  Frequency  Redundancy  Technique  (Upper). 

Resulting  lOOkFIz  Oblique  lonogram  After  Processing  through  the  Modified  Wiener  Trace 
Enhancement  Algorithm  and  the  Noise  Suppression  Technique  (Lower).  Goose  Bay.  Labrador 
to  Millstone  Hill,  Massachusetts:  April  5.  1988,  1908  UT. 


ytm*  :  66  :  096  Tim  :  16:S9:06 


Figure  5.  End  Point  Synthesized  Oblique  lonogram  via  Secant  Law  (Lower).  Region  of  Interest 
Calculated  from  the  E  and  F  Layers  of  the  Corresponding  Vertical  lonogram  (Upper).  Goose 
Bay.  Labrador  to  Millstone  Hill.  Massachusetts;  April  .S.  198K.  18.‘i9  UT. 
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The  sitnulaied  annealing  optimization  technique 
we  use  to  determine  the  QP  electron  density 
profile  has  been  applied  to  solve  many  large 
scale  combinational  optimization  problems  (van 
Laarhoven  and  Aarts,  1988).  In  these  optimiza¬ 
tion  problems  an  objective  function  that  depends 
on  the  configuration  of  a  large  set  of  system  pa¬ 
rameters  is  minimized.  In  the  QP  electron 
density  profile  inversion  problem,  the  system  is 
the  collection  of  quasiparabolic  segments  used  to 
model  the  midpoint  ionosphere.  The  objective 
function  used  in  the  electron  density  inversion  is 
the  sum  of  the  squared  differences  between  the 
group  paths  calculated  from  the  QP  model  and 
those  obtained  from  the  scaled  oblique  iono- 
grams. 

Typical  results  are  seen  in  Figures  6,  7  and  8. 
The  left  panels  of  each  figure  show  the  group 
paths  calculated  from  the  QP  electron  density 
model  (dots)  in  relation  to  the  oblique  ionogram. 
The  right  panels  are  the  corresponding  ground 
ranges  in  relation  to  the  known  ground  range 
(Straight  line).  Seven  scaled  points  from  the 
oblique  ionogram  were  used  in  this  example.  The 
upper  panels  of  Figure  6  contain  the  original 
scaled  oblique  ionogram  and  the  known  ground 
range.  The  middle  panel  of  Figure  6  shows  the 
group  paths  and  ground  ranges  of  the  seven  rays 
for  the  initial  estimates  of  the  layer  parameters 
and  rays  takeoff  angles.  Any  group  path  differ¬ 
ences  greater  than  500km  are  plotted  as  being 
equal  to  SOOkm,  hence  the  four  group  path 
points  at  the  top  of  the  left  middle  panel  of  this 
figure.  After  the  first  iteration  the  three  low 
angle  ray  group  paths  are  within  41km  of  the 
actual  oblique  data,  while  the  three  low  angle 
ray  ground  ranges  are  within  43km  of  the  actual 
ground  range.  Figure  7  contains  the  data 
obtained  for  iterations  3  through  5  in  the  upper, 
middle  and  lower  panels  respectively,  while 
Figure  8  contains  the  corresponding  data  for 
iterations  6  through  8.  As  the  iterations  are 
performed  the  algorithm  converges  to  a  QP  layer 
that  minimizes  the  total  difference  between  the 
measured  and  calculated  oblique  ionograms.  By 
the  eighth  iteration  (lower  panel  of  Figure  8)  the 
largest  group  path  difference  between  the  actual 
data  and  the  quasiparabolic  model  is  4km,  while 
the  largest  ground  range  difference  is  6km. 


4 .  CONCLUSION 

A  small  low  cost  ionospheric  sounder  has  been 
developed  and  tested  that  is  optimally  suited  for 
automated  vertical  and  oblique  sounding. 
Automatic  scaling  algorithms  for  the  vertical  as 
well  as  the  oblique  ionograms  provide  the  verti¬ 
cal  electron  density  profiles  at  the  sounder  loca¬ 
tions  and  at  the  midpoints.  While  the  real  time 


scaling  profile  inversion  of  the  vertical  iono¬ 
grams  has  become  a  routinely  applied  technique, 
the  automatic  scaling  of  oblique  ionograms  in 
terms  of  quasiparabolic  midpoint  profiles  is  still 
at  an  experimental  stage. 
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Figure  6.  Algorithm  Convergence  Properties,  Iterations  0  through  2. 


Figure  7.  Algorithm  Convergence  Properties.  Iicraiions  ?  through  .S. 
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Figure  8.  Algorithm  Convergence  Properties,  Iterations  5  through  1. 
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DISCUSSION 


C.  GOUTELARD 

Voire  conununication  esi  ires  inimssantc  ei  eUe  m'a  passionne  au  poini  de  vous  poser  deux  (juesuons. 

Vous  signalez  que  la  puissance  ^uivalenie  en  impulsion  de  votre  systime  esi  de  8MW.  U  me  semble  cependani  que 
vous  effectuez  un  iraiiemeni  posi-correlaiion  par  iransform^e  de  Founer.  II  faudiait  lenir  compie  d'un  traiiemeni 
similaire  dans  le  cas  du  sysi6me  impulsionnel.  De  plus  en  tdldddieciion  ou  b  fonciion  d'ambiguiid  du  signal  doii  etre 
pnse  en  compie  les  codes  continus  se  trouveni  pdnalisds.  Pouvez-vous  faire  un  commentaire? 

Ma  deuxiime  question  pone  sur  le  modile  que  vous  utilisez  pour  I'ionosphtre.  un  modble  5  irois  couches  quasi 
paraboliques.  Cela  permet  de  metier  des  calculs  analytiques.  Nous  utilisons  un  modile  proche  de  cclui  de  Bradley 
Dudeney  avec  2  couches  quasi  paraboliques  et  1  couche  quasi  liniauc.  Nous  menons  comme  vous  des  calculs 
analytiques  -  par  exemple  avec  des  iquaiions  de  Lagrange  -  pour  modiliser  la  propagation.  Ouel  esi  voire  avis  sur 
I'opponuniti  d'utiliser  I'un  des  modiles  plui6i  que  I'auire? 

your  paper  is  very  interesting  and  it  impressed  me  so  much  that  I  should  like  to  ask  you  tn'o  questtons: 

You  state  that  the  equivalent  pulse  power  of  your  system  is  8MW.  It  seems  to  me  however,  that  you  carry  out  post¬ 
correlation  processing  using  a  Fourier  trantform.  Similar  processing  would  need  to  be  allowed  for  in  the  pulsed 
system.  What  is  more,  the  continuous  codes  are  penalised  in  remote  sensing,  where  allowance  needs  to  be  made  for 
the  ambiguity  function  of  the  signal.  Would  you  care  to  commeni  ? 

My  second  question  concerns  the  model  which  you  use  for  the  ionosphere,  a  model  with  three  quasi  parabolic  layers 
This  enables  one  to  make  analytical  calculations.  We  use  a  model  similar  to  the  one  produced  by  Bradley  Dudeney 
with  2  quasi  parabolic  layers  and  I  quasi  linear  layer.  Like  you,  we  make  analytical  calculations  •  using  Lagrange 
equations  for  example  -  to  model  propagation.  What  do  you  think  is  the  advisabilirv  of  using  one  model  rather  than 
the  other  ? 

AUTHOR'S  REPLY 

1.  For  oblique  sounding  we  use  the  127  bit  maximal  length  sequence  code  providing  21  dB  S/N  improvement.  Another  21 
dB  5/N  gain  is  obtained  by  the  128  point  Doppler  integration.  42  dB  over  500  W  corresponds  to  8  hfW  of  a  single  pulse  no- 
integration  system.  Without  spearal  integration  (like  in  the  Digisonde)  one  should  not  phase-coherently  integrate  more  than 
8  pulses  (to  avoid  loss  of  Doppler  shifted  echoes).  This  requires  1  MW  for  the  same  5/N  radio. 

2.  The  Bradley-Dudeney  profile  may  have  an  advantage  over  the  3  parabola  model  in  that  it  has  less  harsh  discominuities 
between  the  layers.  We  prefer  the  parabolic  model  since  it  can  be  made  completely  smooth  (continuous  first  derivative)  by- 
adding  two  inverted  parabolas  between  the  layer  parabolas.  Fast  analytical  ray  tracing  through  parabolic  sections  of  the 
ionosphere  exists.  Stability  of  the  inversion  solutions  is  good  for  our  technique,  I  do  not  know  what  it  would  be  for  the  B-D 
profile. 


E.  SCHWEiCHER 

1.  Did  you  use  a  Barker  code  or  some  other  pseudo-stochastic  code? 

2.  Did  you  perform  your  pulse  compression  by  software  or  by  hardware? 

AUTHOR’S  REPLY 

1.  A  complementary  8  bit  code  is  used  for  vertical  sounding  in  a  'moderate'  ionosphere.  The  Barker  code  is  used  in  a  high- 
Doppler  (polar  cap)  ionosphere. 

2.  Software  pulse  compression. 
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APPENDIX  B 


Operator's  Manual  for  DORIS  Profile  Inversion 


UNINTERSITY  OF  MASSACHUSETTS  LOW'ELL 
CENTER  FOR  ATMOSPHERIC  RESEARCH 
450  Aiken  Street 
Lowell,  MA  01854 


OPERATOR'S  MANUAL  FOR  DORIS  PROFILE  INVERSION 

April  1992 
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Running  the  DORIS  profile  inversion  software  requires,  as  a  minimum,  the 
following  PC  system: 


Central  Processing  Unit: 
Math  Co-Processor: 
Memory  Size: 

Hard  Disk  with  at  least  4 
Video  Adapter: 

Monitor: 

or  equivalent 
Operating  System: 


Intel  80386  or  higher 
Intel  80387  or  higher 
640  kb  or  larger 
free  Mb 

Color  Graphics  with  VGA  Mode 
NEC  Multisync 

DOS  3.0  or  higher 


The  ULCAR  supplied  ADEP  (ARTIST  Data  Editing  and  Printing)  system 
fulfills  this  requirement. 


To  use  the  off-line  DORIS  SOFTWARE,  do  the  following: 

Go  to  the  DORIS  directory  on  the  'C  drive  by  typing  after  the  prompt: 

C:\  <a> 

CD  DORIS. 


At  this  point,  if  you  type  DIR,  you  will  get  a  listing  of  all  the  files  in  this 
DIRECTORY.  The  executable  files  have  the  extension  .EXE.  There  are  3 
executable  files:  GENSW.EXE,  OBL_INVR.EXE  and  PLOTSYN.EXE. 

The  OBL_INVR  program  produces  mid-point  electron  density  profiles  from 
scaled  oblique  ionograms.  To  run  this  program  type  after  the  prompt: 

OBL.INVR 

then  you  will  be  prompted  to  input  the  scaled  trace  data  file.  A  sample  of 
this  file  is  named  SCLDATA.  Type  SCLDATA. 

The  inversion  algorithm  requires  at  least  2  scaled  points  for  each  mode 
(high  ray  and  low  ray)  of  each  layer  E,  FI,  F2.  resulting  in  12  as  the 
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modes/layer)  (3  layers).  The  first  line  of  this  file  is  the  number  of  scaled 
points.  The  subsequent  lines  are  the  scaled  points.  The  first  column  is  the 
sounding  frequency.  The  second  column  is  the  absolute  group  path.  The 
third  column  is  either: 


1 

for 

low 

ray 

E 

- 1 

for 

high 

ray 

E 

') 

for 

low 

ray 

FI 

-2 

for 

high 

ray 

FI 

3 

for 

low 

ray 

F2 

-3 

for 

high 

ray 

F2 

A  listing  of  this  file  is  enclosed. 

The  program  will  then  prompt  you  for  the  output  file  name.  Type  after 
the  prompt  the  file  name.  Type  (in  this  case  the  output  is  demo  1)  DEMO  1. 

Then  enter  the  ground  range  for  oblique  data.  The  ground  range  for  the 
sample  data  is  2060.00  km. 

Type  2060.00 

Finally,  you  will  be  prompted  to  input  the  data  file  containing  the  mid¬ 
point  quasiparaboiic  electron  density  estimate.  A  sample  file  is  provided 
called  EDPOl.INT. 

Type  edpOl.int 

The  edpOl.int  has  3  columns  and  3  rows.  The  first  column  corresponds  to 
the  E-layer  data.  The  second  column  corresponds  to  the  FI  layer  data.  The 
third  corresponds  to  the  F2  layer  data.  On  the  other  hand  the  first  row  is 
the  critical  frequency.  The  second  row  is  the  semi-thickness  and  the  third 
row  is  the  base  height. 
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Before  the  program  proceeds,  it  will  give  you  a  chance  to  change 
parameters  listed  in  the  electron  density  file  (in  this  case  the  sample  file 
EDPOl). 

Type  0  for  no  change  and  1  to  change  the  parameters.  To  change  the  F2 
layer,  type  1.  After  the  prompt,  type  the  new  parameters.  The  program  is 
now  running.  It  finds  the  updated  quasiparabolic  model  for  the  E  layer.  FI 
layer  and  F2  layer.  The  results  are  dumped  to  the  screen.  Also  8  files  are 
generated.  These  files  are:  DEMOl.edO  and  DEMOl.iOO  containing  the  mid¬ 
point  electron  density  profile  and  the  group  path  of  the  sounding 
frequencies  for  the  initial  input. 

DEMOl.edl  and  DEMOl.iOl  contain  the  mid-point  edp  and  group  paths  at 
the  indicated  sounding  frequencies  after  E  layer  inversion.  DEM02.ed2  and 
DEM01.i02  contain  corresponding  data  after  E  and  FI  layer  inversion, 
DEMOl.edS  and  DEMOl.iOS  contain  the  corresponding  data  after  E,  FI  and 
F2  layer  inversion. 

Another  executable  file  is  gensyn.EXE.  This  program  produces  an  oblique 
ionogram  from  a  three  layer  quasiparabolic  electron  density  profile  file 
specified  by  the  user.  The  ground  range  and  frequency  limits  of  the 
oblique  ionogram  are  also  selectable  by  the  user. 

To  run  this  program,  type  gensyn  after  the  prompt.  You  will  be  prompted 
the  input  the  desired  ground  range.  Type  2060.00,  then  you  will  be  asked 
to  input  the  mid-point  electron  density  data  file.  Type  DEMOl.edS  (note 
that  this  file  contains  the  optimal  mid-point  E,  FI  and  F2  layers). 

The  program  will  give  you  a  chance  to  change  the  input  parameter.  Type  0 
for  no,  1  for  yes.  It  will  also  give  you  a  chance  to  change  the  frequency 
range  data. 

Then  it  will  ask  for  the  output  file  name.  Type  Sample. 

The  third  executable  field  is  plotsyn.EXE.  This  program  plots  both 
synthetic  oblique  ionograms  and  scaled  oblique  ionogram  files. 
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To  run  this  program  type  after  the  prompt  plotsyn. 

It  will  ask  you  to  input  the  power  threshold,  type  0.0.  Then  it  will  prompt 
you  to  input  the  name  of  the  file  to  be  displayed  tj'pe.  Sample  (output  of 
gensyn  program). 

The  program  will  display  the  starting  frequency  and  the  ending  frequency 
and  will  ask  you  if  you  want  to  use  these  values.  Enter  1  for  yes  and  0  for 
no.  Type  1  and  then  it  will  ask  you  if  you  like  to  plot  scaled  points.  Enter 
1  for  yes  and  0  for  no.  Type  1. 

Finally,  the  program  will  ask  you  if  you  would  like  to  plot  inversion  data, 
type  1  for  yes  and  0  for  no.  Type  0. 

Remark:  The  request  for  inversion  data  requires  an  ionogram  inversion 

file  such  as  DEMOl.iO*. 


Then  an  ionogram  is  produced.  The  horizontal  axis  is  the  frequency  axis 
and  the  vertical  axis  is  the  ground  range  axis.  The  ground  range  axis  has 
values  between  2,000  km  and  3.300  km.  The  horizontal  axis  is  shown  in  1 
MHz  steps. 

The  E  layer  is  plotted  in  dark  blue,  1  hop,  and  light  blue  2  hop.  The  FI 
layer  light  and  dark  red  for  1  and  2  hop  respectively.  The  F2  layer  is 
green  with  the  same  light  and  dark  1  and  2  hop  tagging. 
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C:\>d 


l':\>c()  ooTxs 
D:\l>ORIS>d(lu 

SdR  for  DOS  A.e  ind  5,6,  L.  F.  kitrosser,  19S£ 
VoluM  in  drive  D  is  HS-ftOnDRIVL 
Directory  of  D:\D0RIS 


• 

(DIR)  64-14*12 

6:26a  3  PLDTSYH 

EXE 

56682  89-63-91 

2:13p 

•  « 

(DIR)  64-14-92 

6:26a  3  PLOTSYH 

FOR 

5383  69-63-91 

2:13p 

DEIfOl 

36515  64-13-92 

2:52p  3  OPSUBHL 

FOR 

9142  66-26-91 

4;4bp 

DORIS 

DOC 

3256  84-13-92 

16: 16a  3  OPSUBHL 

OBJ 

15366  67-23-91 

4:16p 

EDPei 

IKT 

99  64-89-92 

16:29a  3  ROOT 

EDI 

99  64-18-92 

9:3?a 

GEHEDP 

FOR 

1731  64-69-92 

16:66a  3  ROOT 

ED2 

99  64-16-92 

9:37a 

GEHSYN 

EXE 

45532  84-15-92 

9:36a  3  ROOT 

ED3 

99  64-16-92 

9:37a 

GENSTN 

FOR 

4524  64-13-92 

9:37a  T  ROOT 

106 

526  64-16-92 

9:37a 

6ENSYN 

OBJ 

6476  64-13-92 

9:36a  3  ROOT 

101 

52B  64-16-92 

9;37a 

HOKE 

FOR 

2354  64-67-92 

2:45p  3  ROOT 

102 

564  64-16-92 

9:37a 

NOHE 

OBJ 

3536  64-87-92 

2:45p  3  ROOT 

103 

564  64-16-92 

9:37a 

0BL_INVR  EXE 

56728  64-16-92 

9:47a  3  SCLDATA 

564  67-22-91 

3:24p 

OBL  JHVR  FOR 

13M7  84-16-92 

9:46a  3  SYTEST 

47525  64-13-92  l6:B2a 

OBL  JNVR  OBJ 

16742  64-16-92 

9:47a  3 

Above  display 

324673  bytes  in 

use  27  file(s) 

594944  bytes  free 

D:\DORIS>oblJnvT 
input  sciled  trace  data  file 
scldata 

input  nue  of  output  files 
draol 

input  ground  range  for  oblique  data 

eese.ee 

input  eidpoint  electron  density  estiaate  data  file 
edpBl.int 
present  qp  udel 


E 

fc 

E  rb 

E  y» 

FI 

fc 

FI  rb 

F2  ye 

FE 

fc 

FI  rb 

F2  ye 

4.B6 

96.ee 

26.66 

L.ee 

i2e.ee 

7e.ee 

ie.se 

216.66 

ise.ee 

do 

n 

you  nant 

to  change 

E  layer  yes=l  no=e 

V 

do 

c 

you  «ant 

to  change 

FI  layer  yes-1  no=8 

c 

do 

1 

you  nant 

to  change 

r2  layer  yes^l  no=e 

a 

input  F2  fc  , 

,  FE  H  , 

F2  ya 

ie.:.Eie.e,ise.e 
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Of  It  scilct  points  tne  initiil  estiMt*  Mtcnctf  IS 

Uitn  M  pvpTPff  spttpred  E  froup  path  error  of  112. 

an  avera|e  souared  FI  froup  path  error  of  14S.31S 

an  averafe  squared  F2  froup  path  error  of  34.327 

qp  oodel  updated  durinj  E  la^er  processinf 


E  fc 

Fl 

fc 

F2 

fc 

E  ya 

Fl 

yn 

F2 

y« 

E  rb 

Fl 

rb 

F2 

rb 

4.85 

8.68 

18.38 

IS.  14  7B.M  ite.ee 
S3.57  128.88  218.88 

Of  It  scaled  points  the  current  estioate  oatched  15 

Uith  an  averafe  squared  E  froup  path  error  of  8.882 

an  averafe  squared  FI  froup  path  error  of  177. 48S 

an  averafe  squared  F2  froup  path  error  of  34.484 

qp  lodel  updated  durinf  E  layer  processinf 


E  fc 

Fl 

fc 

F2  fc 

E  yn 

Fl 

yn 

F2  ya 

E  rb 

fl 

rb 

F2  rb 

4.18 

8.88 

18.36 

18.55  78.86  188.88 

37.87  128.88  218.88 

Of  18  scaled  points  the  current  estinate  natched  18 

Uith  an  averafe  squared  E  froup  path  error  of  7.828 

an  averafe  squared  FI  eroup  path  error  of  145.239 

an  averafe  squared  F2  froup  path  error  of  34.548 


qp  aodel  after  E  layer 

processinf 

E  fc 

Fl  fc 

F2  fc 

E  ya 

Fl  y. 

Fe  yn 

E  rb 

Fl  rb 

F2  rb 

4.18 

8.68 

18.38 

18.55 

78.88 

188.06 

37.07 

128.08 

218.06 

qp  nooel  updated  durinf  Fl  layer  processinf 


E  fc 

Fl  fc 

F2  fc 

E  yn 

Fl  ya 

F2  ya 

E  rb 

Fl  rb 

F2  rb 

4.16 

5.94 

16.36 

1£.55 

67.34 

188.66 

97.07 

121.55 

216.66 

Of  It  scaled  points 

the  current  estinate  natched 

With  an  averafe  squared  E  froup  path  error  of  7.828 

an  averafe  squared  Fl  froup  path  error  of  163.374 

an  averafe  squared  F2  group  path  error  of  34.338 

qp  nooel  updated  during  Fl  layer  processinf 


fc 

Fl  fc 

F2  fc 

yn 

Fl  ya 

F2  ya 

rt 

Fl  rb 

F2  rb 

4.16 

S.73 

18.36 

It.SS 

88.36 

188.66 

97.67 

128.31 

218.68 
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Of  It  iCiiti  ^inu  tM  WTTOTt  eStJMt*  MtCtIMl  It 

liitri  *t,  <¥er«f«  &41MTM  t  imp  Mto  erroT  of  T.828 

on  owFTOft  foiurtd  FI  |mp  poth  ottot  of  It.  396 

in  ivoTife  souiTod  F2  imp  pith  ottot  of  37,452 

qp  oodtl  opPitn)  Purini  FI  liytr  prac«$sin| 


E  fc 

FI  fc 

F2  fc 

E  vn 

n  yn 

F2  yn 

1  Tb 

FI  rb 

F2  rb 

4.16 

5.66 

16.36 

16.55 

67.71 

166.66 

97.67 

125.65 

216.66 

Of  16  sciled  points  the  current  estiute  utched 

With  in  iveriie  squired  E  |mp  pith  error  of 

in  iveri|e  squired  FI  |mp  pith  error  of 

in  iveri|e  squired  F2  |mp  pith  error  of 

qp  nodel  updited  durinq  FI  liyer  processini 


E  fc 

FI 

fc 

F2 

fc 

E  yn 

FI 

r* 

F2 

r» 

E  rb 

FI 

rb 

F2 

rb 

4.16 

16.36 

16.55  66.69  166.66 

97.67  127.47  216.66 

Of  16  sciled  points  the  current  estinite  nitched 

With  in  iveriie  squired  E  |roup  pith  error  of 

in  iveri|e  squired  FI  group  pith  error  of 

in  iverige  squired  22  group  pith  error  of 

qp  nodel  updited  during  FI  liyer  processing 


E  fc 

FI 

fc 

F2  fc 

E  yn 

FI 

yn 

F2  yn 

E  rb 

FI 

rb 

F2  rb 

4.16 

5.72 

16.36 

16.55  66.14  166.66 

97.67  136.63  216.66 

Of  16  sciled  points  the  current  estinite  nitched 

Uith  in  iverige  squired  E  group  pith  error  of 

in  iverige  squired  FI  group  pith  error  of 

in  iverige  squired  F2  group  pith  error  of 

qp  nodel  updited  during  FI  liyer  processing 


E  fc 

FI  fc 

F2  fc 

E  yn 

FI  yn 

F2  yn 

E  rb 

FI  rb 

F2  rb 

4.16 

5.76 

16.36 

16.55 

67.51 

166.66 

97.67 

133.71 

2ie.ee 

Of  16  sciled  points 

the  current  estinite  nitched 

Uith  in  iverige  squired  E  group  pith  error  of 

in  iverige  squired  FI  group  pith  error  of 

in  iverige  squired  F2  group  pith  error  of 

qp  nodel  updited  during  FI  liyer  processing 


fc 

FI  fc 

F2  fc 

y» 

FI  yn 

F2  yn 

rb 

FI  rb 

F2  rb 

4.16 

5.73 

16.36 

16.55 

69,23 

166.66 

97.67 

136.97 

216.66 

16 

7.628 

9.761 

36.387 


16 

7.620 

9.271 

39.172 


16 

7.628 

7.554 

36.546 


16 

7.626 

6.916 

36.586 
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Of  U  tciltt  winti  tnc  eurmt  tctiMtt  aatcfttd  18 

Uitfi  «Ti  <vfT4ye  sutured  E  frov^  error  of  7.828 

an  awrrafe  soured  FI  frood  path  error  of  8.91S 

an  averaie  squared  F2  froup  path  error  of  38.46A 

qp  oodel  upoated  danni  FI  lapn  processini 


E  fc 

Fl  fc 

F2  fc 

E  yu 

ri  yi 

F2  yo 

E  rb 

Fl  rb 

F2  rb 

4.18 

5.88 

18.38 

18.  SS 

89.51 

186.N 

97.07 

134.88 

218.81 

Of  18  scaled  points  the  current  estiaate  oatched  18 

Kith  an  averaqe  squared  E  iroup  path  error  of  7.828 

an  averafe  squared  Fl  froup  path  error  of 

an  averape  squared  F2  proup  path  error  of  38.318 

qp  oodel  updated  durinp  Fl  layer  processinp 


E  fc 

Fl  fc 

F2  fc 

E  yn 

Fl  yo 

n  yo 

E  Tb 

Fl  rb 

F2  Ti 

4.18 

5.74 

18.38 

18.55 

88.87 

186.88 

97.87 

133.98 

216.88 

01  16  scaled  points  the  current  estiaate  oatched  18 

With  an  averape  squared  E  proup  path  error  of  7.828 

an  averape  squared  Fl  proup  path  error  of  !).A18 

an  averape  squared  F2  proup  path  error  of  38.8bA 

qp  oodel  after  E  and  Fl  layer  processinp 


£  fc 

Fl 

fc 

F2  fc 

E  yo 

Fl 

y« 

F2  yu 

E  rb 

Fl 

Tb 

F2  Tb 

4.18 

5.74 

18.36 

1&.S5  88.87  188.88 

97.87  133.98  218.88 

qp  oodel  updated  durinp  F2  layer  procesinp 


E  fc 

Fl  fc 

F2  fc 

Eyo 

Fl  yo 

F2  yt 

E  rb 

Fl  Tb 

F2  Tb 

4.18 

5.74 

ie.18 

18.55 

88.87 

188.91 

97.87 

133.98 

289.95 

Of  1£  scaled  points  the  current  estioate  oatched  16 

Mith  an  averape  squared  E  proup  path  error  of  7.628 

an  averarf  squared  Fj  proup  path  error  of  5.419 

an  averape  squared  F2  proup  path  error  of  19.388 


qp  oodel  updated  durinp  r2  layer  procesinp 


E  fc 

Fl  fc 

F2  fc 

E  yo 

Fl  yn 

F2  y. 

• 

E  Tb 

Fl  Tb 

F2  Tb 

4.18 

5.74 

18.89 

18.55 

88,67 

181.75 

97.87 

;3;'.90 

214.98 
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Of  U  Kiltd  points  the  nrrent  estiMte 
Kith  «n  everaie  mured  E  fronp  peth  erroT 
in  ivetiie  squired  FI  qroup  pith  error 
in  iverije  squired  F2  froup  pith  error 
qp  nodel  updited  during  F2  Uyer  procesing 


E  fc 

FI  fc 

FE  fc 

^  r« 

FI  yn 

FE  yi 

E  rh 

FI  rb 

FE  rb 

4.18 

5.74 

18.88 

lt.S5 

tt.t7 

188.83 

97.87 

133.98 

E18.44 

Mtcncd  It 
of  7.t2t 

of 

of  ie.91E 


J 


Of  It  soiled  points  the  current  estinite  nitched 
Kith  in  iverife  squired  E  group  pith  error  of 
in  iverige  squired  FI  group  pith  error  of 
in  iverife  squired  F2  group  pith  error  of 
qp  nodel  updited  during  F2  Uyer  procesing 


It 

7.62B 

5.419 

LL5B 


E  fc 

FI 

fc 

FE  fc 

E  yn 

FI 

y« 

FE  yn 

E  rb 

FI 

rb 

FE  rb 

4.18 

5.74 

18.87 

It.SS  te.t7  17E.26 

97.87  133.98  M9.43 


Of  It  sciled  points  the  current  estinite  nitched 

With  in  iverige  squired  E  group  pith  error  of 

in  iverige  squired  FI  group  pith  error  of 

in  iverige  squired  FE  group  pith  error  of 

qp  nodel  updited  during  FE  liyer  procesing 


E  fc 

FI  fc 

FE  fc 

E  yn 

FI  yn 

FE  yn 

E  rb 

FI  rb 

FE  rb 

4.18 

5.74 

18.13 

ie.S5 

8t.67 

189.58 

97.87 

133.98 

E13.5B 

Of  It  sciled  points  the  current  estinite  nitched 

With  in  iverige  squired  E  group  pith  error  of 

in  iverige  squired  FI  group  pitt  error  of 

in  iverige  squired  FE  group  path  error  of 


qp  nodel  ifter  E  FI  ind  FE  liyer  processing 


FI  fc 

FE  fc 

I 

FI  yn 

F2  y» 

1 

FI  rb 

FE  rb 

4.18 

5.74 

18.13 

18.55 

8t.t7 

189.58 

97.87 

133.98 

E13.SB 

It 

7.tse 

5.419 

7.60E 


IB 

7.BEe 

5.419 

t.t74 
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(■:\CORlS^f«iKm 
in»ttt  tesirti!  iroum)  r<n|t 

zoe.i 

iR^  ixofoiiit  clertm  wmity  iiU  fiit 
dtMl.c43 


present  qp  oodel 

E  fc  E  rb 

E  VO 

FI  fc 

FI  rb 

F2  yo 

F2  fc 

FI  rb 

F2  ye 

4.18 

97.87 

16.55 

5.74 

133.98 

66.67 

18.13 

213.58 

169.58 

do  you  Mot  to  chinfc  E  liyor  yM»l  oo^ 

B 

do  you  Miot  to  chonfe  FI  lnyoT  yos*!  i»*8 

8 

do  you  Mot  to  chin|e  F2  Uyer  yes=l  i»«8 

8 

pretent  freduency  T*n|e  dita  fo^dltfifoai 

:;.8e8  .lee  28.888 

do  you  Mint  to  chanie  ttiete  vilues  yes*l  bo*8 

8 

leneretiny  ofelique  lonoftM 

5.888880888888898  1.8888888 1^981 lElE-881 

input  Mae  of  output  file 
tteplc 


0 


l:\D0RIS) 
l':\DORIS)plotsyii 
input  pouer  tnreshold 


18.138888 
189.588888 
213.5&B888 

fein  »  5.8B9e8efei»  *  28.888888 

do  you  uent  to  use  these  vslues  l=yes  B^no 
1 

do  you  uint  to  plot  scaled  points  I'yes  8-no 
1 

inout  nsee  of  scaled  point  file 
scloata 

oc  you  uant  te  plot  inversion  data  l*yes  8=no 
8 


V 

input  naoe  of  file  to  he  displayed 
saeplq 

4.188888  5.748888 

It. 558888  68.678888 

97.878888  133.988888 
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